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Abstract 

Superatomic molecular orbitals (SAMO) in Cgo are ideal building blocks for functional nanos¬ 
tructures. However, imaging them spatially in the gas phase has been unsuccessful. It is found 
experimentally that if Cgo is excited by an 800-nm laser, the photoelectron casts an anisotropic 
velocity image, but the image becomes isotropic if excited at a 400-nm wavelength. This diffuse 
image difference has been attributed to electron thermal ionization, but more recent experiments 
(800 nm) reveal a clear non-diffuse image superimposed on the diffuse image, whose origin remains 
a mystery. Here we show that the non-diffuse anisotropic image is the precursor of the / SAMO. 
We predict that four 800-nm photons can directly access the If SAMO, and with one more photon, 
can image the orbital, with the photoelectron angular distribution having two maxima at 0° and 
180° and two humps separated by 56.5°. Since two 400-nm photons only resonantly excite the 
spherical Is SAMO and four 800-nm photon excite the anisotropic If SAMO, our finding gives a 
natural explanation of the non-diffuse image difference, complementing the thermal scenario. 
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I. INTRODUCTION 


In nanoscience, no other cluster has ever garnered more attention than C 


60 


Q. 


Such 


a highly__sy;mmetric molecule allows a high electjron^ delocalization, enabling ultrafast c^- 


J, l5| , and high harmonic generation la la. 

discovered a 


namics { 2 , 3|, strong nonlinear optical responses 

Recently, employing scanning tunneling microscopy (STM), Feng et al. 
group of unusually large orbitals both inside and outside of the Cgo cage, which are both 
beautiful and surprising. These orbitals, called superatomic molecular orbitals (SAMOs), 


have attracted immediate attention 


10| . Spatially, they bear a close resemblance to their 


atomic counterparts, but with a much larger radius. This motivated Roy and coworkers 111 
to design nanoscale atoms for solid-state chemistry. Figure [U shows a 1/ orbital, where one 
sees a distinctive anisotropy characteristic of / orbitals. However, imaging such high-lying 
orbitals becomes increasingly difficult for STM, because of the strong overlap between the 
orbitals of Ceo and those of the substrates. 

The velocity-map imaging technique 


12 


13| does not have this problem, since it works in 


the gas phase. Figure [U schematically shows that after the laser strikes Cgo, the photoelec¬ 
trons with different velocities, after several stages of accelerating electrode plates, cast an 
image on the phosphor screen, which is captured by the CCD camera. This image carries the 
brmation. Doing so, quite surprisingly, neither Johansson et al. nor Kjellberg 
detected any image similar to superorbitals. Instead, they reported a seemingly 


orbital in 
et al. \l4 


irrelevant observation: The electron image is anisotropic if excited by an 800-nm laser, but 
is isotropic with a 400-nm laser. How could it be possible that a similar observation in ion 
yields (ellipticity dependence) [3, 1^ occurs for the electrons as well? Huismans et al. 171 
further showed that the anisotropy increases with the laser intensity, and they assigned it 
to the accumulation of highly structured but slightly different angular distributions of the 
ground states, while others assign the diffuse part of the photoelectron spectrum to thermal 
electron ionization or a “kick” from the laser field lOj, with an explanation presented in 
the second paper of Ref. [lOj. Huismans et al. emphasized that a fully satis; 
requires a complete ab initio calculation. Such a calculation is now available 
studies did not address the difference between the 800-nm and 400-nm excitations. Very 
recently, Li and coworkers experimentally showed a distinct and non-diffuse image superim¬ 
posed on the diffuse anisotropic background |^. These non-diffuse features can not easily 


actory theory 
1^ . but these 
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be explained by the thermal electron ionization, since the Boltzmann-like distribution is 
unlikely to yield an image with intensities only concentrated at six spatial locations. We 
wonder how such nondiffuse images are formed, besides the origin in the thermal electron 
ionization. 

In this paper, we demonstrate theoretically that SAMOs in Cgo in general, and the If 
orbital in particular, are accessible to the velocity map imaging technique. Our first im¬ 
portant prediction is that the energy gap between the 1/ orbital and HOMO matches the 
four-photon energy of an 800-nm laser pulse. By successively absorbing four photons, elec¬ 
trons in the HOMO can be excited into the 1/ orbital; by absorbing one extra photon, they 
are cleared of the ionization potential energy of 7.58 eV, and can then be detected by the 
velocity map imaging. We further predict that the angular resolved spectrum has a distinc¬ 
tive If orbital feature, with two local maxima, separated by 56.5°. There are two reasons 
why the If orbital can be effectively probed. First, the multiphoton excitation suppresses 
an otherwise strong contribution from nonSAMOs. Second, although energetically d and g 
SAMOS are in the vicinity of the 1/ orbital, due to the selection rule, transitions to those 
orbitals are not possible. Such a detection scheme is quite generic. By slightly tuning the 
photon energy of the laser beam, it is also possible to observe 2p, 3p and Ih orbitals. If 
Ceo is excited by two 400-nm photons, the strongest transition is between HOMO-1 and 
nonSAMOs, but this transition is off-resonant by 0.62 eV. On the other hand, although the 
transition matrix element product is small for the transition from HOMO-2 to Is SAMO, 
this transition is nearly resonant, and it is likely that the Is SAMO is excited, consistent 


0 - 


with Ref. lldl. Therefore, our study suggests a new and alternative explanation to the 
anisotropic and non-diffuse image 19|], and demonstrates the new power of multiphoton 
emission for investigating SAMOs in fullerenes. 


II. THEORETICAL FORMALISM 


We employ two complementary first-principles methods within the density functional 
theory: the basis function-based methods (GaussianOQ 20| and VASP 2l|]) and the real 
grid mesh method (Octopus) [2^. To build a case to directly probe SAMOs in Ceo, we first 
investigate the density functionals’ effect on the energy gap between the HOMO and If 
SAMO for a group of six different functionals, within the same basis 6-311-|-G(d). Figure 
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Ha) shows that the functional has a substantial impact on the energy gap AE (see the filled 
red circles). Since the ionization potential (IP) is at 7.58 eV (the dashed green line) and 
If SAMO must be a few eV below IP, some popular functionals such as B3LYP and CAM- 
B3LYP predict the results too large even in the single-particle limit. This result is rigorous 
since it does not involve the complicated photoexcitation and ionization. This convinces us 
that the SWVN/LDA functional is reliable, at least within the single-particle limit. 

Next, we employ the same SWVN/LDA functional but increase the basis function size by 
adding more diffuse functions {q] to the Gaussian basis. When we use the aug-cc-pvtz basis, 
we find that AE drops to 7.136 eV (see the empty red circle in Fig. H^))- We then put a 
layer of 60 hctitious hydrogen atoms at a radius of 9 A from the center of Cgo- These atoms 
have no charge and only serve as the centers for additional basis functions. The carbon 
atoms have the 6-311-l-G(d) basis, while those fictitious hydrogen atoms have aug-cc-pV6Z, 
which contains 7s, 6p, 5d, 4/, 3g and 2h Gaussian primitive functions. Doing so, we find 
that the LDA gap AEif^iDA decreases from 7.506 eV (without hctitious atoms) to 6.957 eV 
(with hctitious atoms), with the net change as large as 0.549 eV (compare the hrst two hlled 
circles of Fig. H^))- is convincing that the basis functions strongly ahect the accuracy 
of the energy levels of SAMOs. However, if we continue to use the Gaussian primitives, the 
calculation becomes increasingly demanding. We must pursue a diherent approach. 

To gradually eliminate the basis function ehect, we employ a grid mesh-based method as 
implemented in Octopus 22|]. An important advantage is that the real grid is a balanced 
approach and treats SAMOs and nonSAMOs on an equal footing. Figure shows that as 
r increases (the bottom horizontal axis represents r), AE gradually converges to 6.337 eV 
(see the squares at the bottom of the hgure). We then ht the data to AE = AEq + AEi/r 
and find AEq = 6.015 eV and AEi = 10.0707 eV/A. AEq = 6.015 eV is considered as a 
lower limit for this gap. As an independent check, we employ a planewave basis function 
as implemented in the VASP code (2l|, where we place Cgo in a big fee supercell with two 
respective lattice constants of a = 21 A and 25 A, and two planewave cutoffs of 500 and 600 
eV. The results are shown in Fig. ^^e plus signs). The VASP result is 6.36 eV, 

which is considered as an upper limit. If we average these two limits, we expect that this 
gap settles down at 6.2 eV, which precisely matches the four-photon energy of a laser pulse 
of wavelength 800 nm. Figure Hb) shows the energy spectrum of SAMOs referenced to the 
HOMO. For an 800-nm laser, five photons are needed to reach the IP (see the thick dashed 
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green line on the top). The fourth photon has access to the SAMOs which are congested in 
a narrow window. Figure [2])c) zooms in on the SAMOs. Consistent with Feng’s VASP result 
8|, the gap between the s SAMO and LUMO is 3.25 eV, but our interest is in higher SAMOs. 
This energy gap becomes an outstanding case to judge the quality of the two most popular 
density functionals, LDA versus B3LYP. If the PDA result is correct, then an 800-nm laser 
can detect 1/ SAMO experimentally. This is the hrst testable case for experiments. 


III. TWO-PHOTON VERSUS FOUR-PHOTON EXCITATION 


Regardless of the experimental finding, such a multiphoton excitation requires a stronger 
laser. Would it be possible to do the same job with two UV photons as Johansson et al. [18 1 
did? In order to have a successive optical transition from the occupied states to unoccupied 
states, the accumulative product of matrix elements of each transition must be different 
from zero. Different from prior studies, we compute the transition matrix elements between 
the HOMOs, SAMOs and nonSAMOs by integrating 


/ OO 

dT'il)i{r)r'ijjf{r), ( 1 ) 

-CO 

where "01 (r) is the Kohn-Sham wavefunction, and r is the electron coordinate. For the 
frequency-doubled case Q. at ,eaat two photoas a.e needed to reaeK the 1/ SAMO. Thus, we 
compute the product of two transition matrix elements, (occ|D|i)(i|D|j), for all the relevant 
transitions around this two-photon energy. The results are surprising. Although there are 
lots of states energetically accessible, only a very few actually contribute. Figure [3](a) shows 
that for two-400-nm-photon photoexcitation, two groups of nonSAMOs have the largest 
product of the transition matrix elements, followed by /, s, d, g and h SAMOs. Although 
the nonSAMOs dominate the product of the matrix elements, their energies are off-resonant 
from two-400 nm photoexcitation. Specifically, the energy for two 400-nm photons is 6.2 eV, 
but the excitation from HOMO-1 to the nonSAMOs is 6.83 eV, off by 0.63 eV, which greatly 
reduces the excitation probability into the nonSAMOs. But for the Is SAMO, although its 
matrix product is small, the excitation energy from HOMO-2 to Is SAMO is 6.17 eV, nearly 
resonant with the two-400 nm photon energy of 6.2 eV. Therefore, during 400-nm excitation, 
the Is SAMO is excited strongly. This hnding is also consistent with Johansson’s assignment 
of their first peak to Is SAMO |^. Since the Is SAMO is spherical, this explains why at 
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400-nm excitation the photoelectron image is isotropic, besides the contribution due to the 
diffuse mechanism as discussed above. 


In the same spirit, we investigate the four-photon process. Four-photon excitations 
build upon four successive transitions, with the product of four transition moments as 
{^|D|j)0'|D|fc)(fc|D|0{/|D| m), where |i) is the HOMO, \m) are the SAMOs, and the others 
are intermediate states. Figure [3]^b) plots the products for all the states with large contri¬ 
butions. Our results are insightful. First, the transition matrix element strictly obeys the 
dipole selection rule and has the correct parity. We see that only SAMOs with odd angular 
momentum quantum numbers have nonzero products. For instance, s and d SAMOs have no 
transition. Second, the SAMO dominates the four-photon transitions. The largest product 
is from the Ih SAMO, and the second largest is from the 1/ SAMO, followed by the 2p and 
3p SAMOs. It is likely that they compete with the If SAMO, but which one dominates 
depends sensitively on the actual photon energy, a finding that allows one to probe various 
SAMOs by tuning the laser wavelength. We also notice that the nonSAMO contribution is 
still sizable. Therefore, to observe SAMOs, future experiments must steer the laser wave¬ 
length away from those nonSAMOs energetically and use an optimal laser intensity (a laser 
too strong would populate other nonSAMOs). For an 800-nm laser, 1/ is strongly excited, 
both in terms of the transition matrix elements and the transition energy. Since its shape is 
highly anisotropic (see the far right schematic in Fig. [T]), this explains why experimentally 
Li and coworkers 19| find that the photoelectron image is non-diffusive and anisotro pic, and 
more importantly has fine structures superimposed on the diffuse image seen before lOl. Il4|. 
To test our results quantitatively, we compute the kinetic energy of the electron. To eject 
this electron from If SAMO, one more photon is needed to clear the ionization potential 
energy IP = 7.58 eV. This leaves the ionized electron with a kinetic energy of 5hu — IP = 0.2 
eV. Sure enough. Fig. 2(e) of Ref, 


■ Q 


indeed shows a peak at 0.2 eV, but this important 
finding did not catch the attention of the investigators, as their emphasis was on the 400-nm 
results. 


IV. DISCUSSION 

To really image the If orbitals, we have to overcome two challenges. First, since the If 
SAMO is seven-fold degenerate, would the final image of If SAMO be completely washed 
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out because of the spatial average during the experiment, in particular as Ceo spins rather 
rapidly? This is an issue for the fluorescence spectrum in CO 2 , where the molecule has no 
preferred orientation. It has to be properly aligned before it can be imaged 23(]. Fortu¬ 
nately, we find that the transition paths for seven If SAMOs are different, so their products 
of transition matrix elements are not the same. Thus, no spatial average is necessary. Nu¬ 
merically we find that they range from 3.2 to 2.1 A'^ (see the circles at the / SAMO in 
Fig. |3](b)). To understand the reason behind this, as an example, we take two 1/ or¬ 
bitals with spherical harmonics F 30 and Fsi, respectively. If the light is linearly polarized, 
for 1^30, the transition path can be like |HOMO) —)■ Fqo Fio Y 20 —)■ I 30 , but for 1 ^ 31 , 
|HOMO) |i) —)■ \j) —Yu Y 21 —?• Fsi, where \i{j)) are two different intermediate states. 
Therefore, the experimental geometry (Fig. [T]) intrinsically differentiates among different If 


states, already found in atoms 


24| . In Ceo we find that the dominant path or doorway state 


25l is |HOMO) ^ \ld) \2p) \2d) \lf). 


The second challenge is whether If SAMO can ever cast its distinctive image on the hnal 
continuum state. We decide to compute the angular distribution of the electron density for 
an / orbital with the spherical harmonic Y 30 . The intensity is proportional to the square of 


(V'(r)|^e *^'0 = 47 ri^h/cy 3 o( 6 'fc, 0 fc) / drRf{r)j^{kr)rf 


( 2 ) 


where 9k and fk define the direction of the wavevector k of the ejected electron, and a 
constant term due to the orthogonality is not included. The angular distribution comes 
from YfQ{6k, fk)- Employing the realistic wavefunction of If SAMO in Ceo, we compute the 
angular distribution of the If SAMO. The results are shown in Fig. [3](c). The maxima are 
at 0 ° and 180°, along the positive and negative directions on the 2 ;-axis, respectively (see 
the / SAMO in Fig. [T]). As the angle rotates away from the 2 :-axis, the intensity drops. We 
believe that this produces the anisotropy seen in the experiment 19|. The If SAMO has two 
distinctive double humps, separated by 56.5°. The unpublished experimental photoelectron 


angular distribution by Li and coworkers 


19j is shown in Fig. [3](d). The agreement between 


our theory and their experiment is remarkable, where the shape, peak and hump locations 
are identical. We draw two dashed lines to highlight the locations of two humps; and our 
estimated experimental angle separation between them is 54.4°, in quantitative agreement 
with our theory within the experimental error. Figure 4.6(f) of Ref. 19|] looks more like our 
predicted If SAMO (see Fig. [ 1 ]). This constitutes a strong experimental support for our 
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theory. 


V. CONCLUSION 


While much attention has been given to the diffuse spectrum in VMI, recent experimen¬ 
tal evidence shows a signihcant contribution from the non-diffuse photoelectrons 19|]. Here 
we demonstrate through a series of carefully designed first-principles calculations that the 
non-diffuse and anisotropic velocity image observed at 800 nm is likely due to the excitation 
from the highest occupied molecular orbital into the If SAMO. We predict that the energy 
gap between the highest occupied molecular orbital and 1/ SAMO matches the four-photon 
energy of 800 nm, and a peak at 0.2 eV will show up in the photoelectron spectrum. In 
agreement with the prior study, two 400-nm photons can excite Is SAMO from HOMO-2 
nearly resonantly, although nonSAMOs have a larger transition matrix product, but they 
are off resonance by 0.6 eV. Since the Is SAMO is spherical, the photoelectron image is 
isotropic, in addition to the possible thermal and diffuse electron emission. Future experi¬ 
mental verification will have a series of important theoretical consequences. Which density 
functional, LDA or GGA, is more accurate to describe the photoelectron angular distribu¬ 
tion and multiphoton emission? Which basis functions, Gaussian diffuse functions or real 
grid mesh, are more accurate to describe superorbitals? For the first time, we demonstrate 
that it is the multiphoton excitation that greatly boosts the transition probability into these 
superorbitals, a hnding that is expected to inspire new experimental and theoretical inves¬ 
tigations. 
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FIG. 1: Proposed experimental setup of velocity map imaging. Photoelectrons from Ceo pass 
through electrode plates and form images on a phosphor screen, which are in turn captured by a 
CCD camera. Depending on whether an 800-nm or 400-nm laser is used, the image is anisotropic 
or isotropic. The diffuse origin of this image difference is attributed to the thermal photoelectron 
emission, but this can not explain the non-diffuse features seen in a recent experiment {l^ . Far 
right; For convenience of viewing, we rotate the / superatomic molecular orbital (SAMO) by 90°. 
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FIG. 2: (a) Energy gap between the If SAMO and HOMO computed with different functionals 

and basis functions (top axis) and as a function of sphere radius (bottom axis). The filled circles 
represent the results computed with Gaussian09, while the empty squares are those obtained with 
Octopus (the pseudopotential and on the real grid mesh). We expect that in the infinite sphere 
limit, the gap settles down at 6.2 eV. The VASP results (plus sign) are included for comparison, 
where the radius refers to the center-to-center distance of Cgo- (b) Energy spectrum of SAMOs 
referenced to the HOMO energy. The ionization potential energy is at 7.58 eV (see the thick dashed 
green line at the top), (c) Zoomed-in SAMO energy spectrum. The If SAMO is between 2d and 
Ig. The SAMOs form the doorway state for the ionization. 
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FIG. 3: (a) Product of two transition matrix elements for the two-photon process at the 400- 

nm wavelength. States that strongly contribute to the signal are highlighted with red ellipses. 
Although the nonSAMOs have a large transition matrix element, due to the off-resonance, their 
contribution is small. By contrast. Is SAMO excitation is nearly resonant, so it dominates the 
photoelectron spectrnm, consistent with the experimental findings, (b) Products of four transition 
matrix elements for the fonr-photon process excited at 800 nm. In contrast to (a), four SAMOs 
(2p, 3p, If and Ih) and one nonSAMO contribute strongly. Since only a subset of SAMOs are 
excited, this casts an anisotropic image, (c) Photoelectron angnlar distribution as a function of 
angle at 800 nm. Two hnmps are separated by 56.5°. (d) Experimental results from Ref. [l^ . 
Excellent agreement with our theory is found. 
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